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ABSTRACT 


At thermionic-diode temperatures, work functions of tungsten electrode surfaces 
depend strongly on positive-ion sheaths between these electrodes and cesium plasmas. 
The present theoretic study demonstrates this effect for a plane field-modified emitter 

with a collisionless positive- ion sheath in a thermally ionized plasma. The analysis 

12 

covers temperatures from 1700 to 2600 K; plasma electron number densities of 10 , 

13 14 3 

10 , and 10 cm” ; and a bare work function of 4. 59 V. Values for sheath voltages, 

work functions, and cesium arrival rates and coverages on tungsten result from equili- 
brium calculations. In addition to these variables, computations for elevated plasma 
electron temperatures yield net current densities. Throughout this work, the sheath 
drops lie below 1 V. And atom arrival rates approach total cesium arrival rates only at 
low temperatures and sheath voltages. 


STAR Category 25 



SOME INTERACTIONS OF POLYCRYSTALLINE TUNGSTEN 
ELECTRODES WITH CESIUM PLASMAS 
by James F. Morris 
Lewis Research Center 

SUMMARY 

Positive-ion sheaths greatly influence surface properties of polycrystalline tungsten 
at thermionic-diode temperatures in cesium plasmas. Indications of this effect on work 
functions come from the present analysis; here, plane field-modified emitters of elec- 
trons, ions, and atoms obtrude collisionless positive-ion sheaths into thermally ionized 
plasmas. For calculations based on these components, temperatures range from 1700 to 
2600 K; cesium plasmas contain 10^, 10^, and 10^ electrons per cubic centimeter; 
and the bare polycrystalline tungsten has a 4. 59- volt work function. Equilibrium sheath 
solutions yield overall voltages, cesium arrival rates and coverages on the electrodes, 
and work functions. These variables as well as net current densities result from compu- 
tations for plasma electrons up to 200 K hotter than the other constituents. 

The generalized sheath properties follow correlations given in previous publications 
on this model: They apply only for equilibria, near equilibria, and relatively low net 
transport. They are functions of the plasma electron number density N g p, the absolute 
temperature of the electrode T^, the work function minus the plasma potential, and the 

IT 1 /r / a\ll V 2 

emission Debye length, A DE (cm) « 6.9^ T E (K) / N ePv m / | * And they result 

from sheaths with effective widths of approximately 2 to 2^ emission Debye lengths. 

Previous studies with parametric work functions revealed no restriction of sheath 
voltages. But the present analysis with work functions affected by cesium adsorption 
produces sheath drops below 1 volt only. This limit derives from restricted ionization 
at the electrodes. Inhibited ion production occurs because of reduced work functions at 
low temperatures and because of decreased cesium arrivals at high temperatures even 
with increased work functions., When the sheath drops as well as the temperatures are 
low, atoms reach the electrodes at rates approaching those for total cesium arrivals. 
Otherwise adsorption of the ions generally outruns that of the atoms by over an order of 
magnitude. Positive- ion sheaths in cesium plasmas definitely react to and strongly 
affect electrode work functions. 



INTRODUCTION 


References 1 to 4 develop theories for plane collisionless sheaths separating ther- 
mally ionized plasmas (ref. 5, Saha) from field-modified emitters of electrons, ions, 
and atoms (refs. 6 to 9: Richardson, Dushman; Saha, Langmuir; and Schottky). These 
models for positive-ion and electron sheaths apply for equilibria, near equilibria, and 
relatively low net transport. The previous investigations (refs. 1 to 3) for parametric 
variations of work functions, temperatures, and plasma electron concentrations give 
sheath characteristics for cesium plasmas. As references 1 to 3 note, however, cesium 
changes electrode work functions and makes coincidence of the parametric and actual 
work functions fortuitous. The present study overcomes this theoretic deficiency. It 
reveals effects of adsorption on work functions of polycrystalline tungsten surfaces with 
positive- ion sheaths in cesium plasmas. 

Collisionless sheaths pass atoms freely but alter ion transport, hence adsorption, 
and therefore surface potentials like work functions. In electron sheaths, emitted ions 
accelerate into the plasma, while ions moving toward the electrode decelerate and often 
reflect. Then the electron sheath reduces the total arrival rate of the plasma chemical, 
ions and atoms (ref. 10). For plasmas with low degrees of ionization, though, electron 
sheaths affect adsorption only slightly. If conditions require solutions for electron 
sheaths, a direct procedure yields them. The sum of atom and ion fluxes leaving the 
plasma serves as an approximate arrival rate for an estimate of the work function 
(ref. 11). With this work function, the emitter temperature, and the plasma conditions; 
references 1 to 3 allow a good guess at the overall sheath voltage. The sheath drop in 
turn yields a correction for the ion arrival rate, then for the work function, and again 
for the sheath drop. Finally a close approximation of the emitter, sheath structure 
evolves. Because of this simple iterative method and the overwhelming interest in posi- 
tive-ion sheaths, the present analysis bypasses the solution for electron sheaths. 

But the positive-ion sheath complicates determinations of arrival rates at electrodes 
(ref. 10). And because this sheath traps many ions at the electrode surface, atom fluxes 
alone from plasmas with low ionization often approximate total arrival rates poorly. In 
positive-ion sheaths, plasma ions accelerate to the electrode, while emitted ions decel- 
erate toward the plasma and often return to the emitter. Such movements oppose those 
in the electron sheath, where total arrival rates depend only on the sheath drop and the 
ion and atom fluxes from the plasma. In addition to their dependence on these three var- 
iables, however, total arrival rates for positive-ion sheaths are functions of three more: 
the ionization potential of the plasma chemical; the sheath field at the electrode; and the 
work function, which varies with surface coverage (ref. 12). The present method com- 
bines the total arrival rate with a work function correlation (ref. 11) in another iterative 
loop on the original solution for positive-ion sheaths (refs. 1 to 4). Now the machine 
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solves for the cesium-ion sheath properties for a specified plasma and bare polycrystal- 
line tungsten, obtains the total cesium arrival rate for that solution, corrects the work 
function for the cesium coverage on the tungsten electrode, computes the sheath proper- 
ties for the new work function, and continues to iterate to a final answer. Without output 
plots (refs. 1 to 4), a complete solution requires about 0. 03 minute on an IBM 7090. 

In the previous discussion relating cesium arrival rates to work functions, refer- 
ences 11 and 12 appeared. Many other papers give theories for adsorption effects on 
surface potentials, but these two are unique. Reference 12 presents perhaps the best 
general solution for this problem thus far; it relates surface properties to coverages. 

As reference 10 indicates, though, desorption times and energies for ions and atoms are 
prerequisites for determinations of surface coverages from arrival rates. These de- 
sorption data vary with coverage and result only from well- controlled experiments. At 
present nothing approaching a complete set of desorption characteristics for an adsorbate 
on a particular crystal face exists. And although Taylor and Langmuir (ref. 13) pub- 
lished a monumental study of cesium on polycrystalline tungsten long ago, technology now 
demands definition of crystal orientation. In the absence of better information, however, 
this report analyzes cesium and polycrystalline tungsten. Furthermore, since refer- 
ence 11 correlates three sets of cesium and tungsten data admirably (refs. 13 to 15), 
even though its physical model is questionable, this work uses that correlation for its 
work function changes. Reference 12 also handles Langmuir’s data well but not nearly 
as precisely or as easily as reference 11. 

Further clarification of the terminology and contribution of reference 11 seems nec- 
essary at this point. First, ’’polycrystalline" implies a surface composed of patches of 
various crystal faces. If these patches are very small compared with the sheath width, 
the surface is acceptable. This requirement provides the essentially uniform surface 
potential assumed in the present sheath model. Second, the assumption of a 4. 59-volt 
work function correlates information from references 13 to 15. But the work functions 
for those references undoubtedly differ from 4. 59 volts: Langmuir used a tungsten wire 
with predominantly 110 surfaces, Houston studied a polycrystalline tungsten, and 
Breitwieser worked with nearly monocrystalline tungsten having a work function of 
5.0 volts. In any event, reference 11 correlates the data of both Langmuir and Houston, 
both of whom reported work functions very near 4. 59 volts. 

Thus the present model describes a plane positive- ion sheath separating a thermally 

ionized cesium plasma from polycrystalline tungsten electrode. The work function 

varies with surface coverages of the emitter below its bare value of 4. 59 volts. Of 

course cesium adsorption on the tungsten depends on the system properties selected for 

the analysis: The temperatures rim from 1700 to 2600 K with plasma electron tempera- 

12 

tures up to 200 K higher than the others; and the plasma concentrations result from 10 , 

13 14 

10 , and 10 electrons per cubic centimeter. From these parameters many detailed 
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results emerge (refs. 1 to 4), but only those that show the effects of the variable work 
function appear here. Sheath voltages, cesium arrival rates and coverages on tungsten, 
and work functions alone indicate equilibria. Added to these variables, net current den- 
sities reveal the near-equilibrium influences of elevated temperatures for plasma elec- 
trons. 

Because the basic sheath model remains unchanged (refs. 1 to 3), the previously 
published correlations of general sheath properties apply to the present study: Dimen- 
sionless sheath widths, voltages, and emitter fields and densities of the majority charge 
at the electrodes depend on the number of electrons in a cubic centimeter of plasma N g p, 
the kelvin temperature of the emitter Tg, and the work function minus the plasma poten- 
tial. Since emitted particles dominate this sheath, the plasma Debye length yields as a 
correlating parameter to the emission Debye length (refs. 1 to 3): 


A DE 



( 1 ) 


Most of the sheaths have widths between 2 and 2-^ emission Debye lengths at these condi- 
tions. In that distance, though, the ion concentration can increase by more than a factor 
of 100 across a 1-volt sheath. Thus fractional ionization can be low in the plasma but 
high next to the emitter. Again the positive- ion sheath stands out as a prime influence 
on cesium arrival rates at electrodes. 


SYMBOLS 

A,B constants in equation for cesium vapor pressure 
E sheath field 

e electronic charge (in dimensions consistent with equations in which it is used) 

I ionization potential 

j current density of particle type 

m particle mass 

N number density of particle 

p vapor pressure 

T absolute temperature 

AV voltage above plasma potential 
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9 fractional surface coverage 

k Boltzmann constant 

X shielding length 

ji particle arrival rate 

A cp change in work function 

<p work function (without subscript); plasma potential relative to electrode Fermi 
level (with subscript) 

Subscripts: 

a atom or atomic 

Cs cesium 

D Debye 

E emission or at emitter 

e electron or electronic 

i ion or ionic 

o plasma potential subscript for near- equilibrium positive-ion sheaths 

oo plasma potential subscript for equilibrium sheaths 

P in plasma 

S across sheath 


THEORY 


References 1 to 4 detail the thoery for plane collisionless sheaths between emitting 
electrodes and thermally ionized plasmas. The electrodes emit electrons, ions, and 
atoms that follow the Richardson, Dushman; Saha, Langmuir; and Schottky relations. 
And the plasma complies with the Saha equation. References 1 to 3 discuss limitations 
of the sheath model and generalize results for cesium plasmas and parametric work 
functions. The present work incorporates a polycrystalline tungsten emitter with work 
functions that depend on the amount of its surface covered by cesium. 

For the positive-ion sheath, the arrival rate of the plasma chemical is 


^aP + ^iP + hE 


M = 


eAV, 


1 - exp 


S> 


k T 


E, 


( 2 ) 
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where AVg is a positive number. Ionic and atomic current densities (j 
pass unchanged from the plasma through this sheath to the emitter. Because no expo- 
nential cutoff of the ion flow occurs, equation (2) imposes no particular form of velocity 
distribution on the ionic and atomic streams from the plasma. But current densities of 
emitted ions j^ E experience exponential cutoffs and particularity reflect back to the 
emitter surface. Therefore, equation (2) requires half- Maxwellian ionic emission, 
which is what the thermal emitter provides. Thus, equation (2) stands essentially unre- 
stricted for cesium on tungsten. 

For the present model of a positive-ion sheath, ionic emission involves the ionization 
potential I and current densities of the plasma atoms and ions in addition to the voltage 
across the sheath and the electrode properties: the absolute temperature T E , work 
function <p, and sheath field E E (refs. 1 to 3): 


iP and j a p) 


*iE “ 


^aP + •'iP 


2 exp 


e(I - <P) + (e 3 E £ ) 
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K T 
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eAV c 
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(3) 


Substituting this expression into equation (2) yields equation (9b) of reference 10: 


M = 


j aP + j iP 
e 


1+2 exp 

/ 3 N 1 / 2 

e(I - <P) + (e^Egj 


« t e 

( eAV s\ 

pvn | 4. 

2 exp 

/ 3 X 1 / 2 

e(I - <p) + (e E E j 

V “V 

kT e 


(4) 


Total arrival rates come from equation (2) because of the use of results from the 
calculations for positive-ion sheaths (refs. 1 to 4). Both ionic and atomic cesium ar- 
riving at the polycrystalline tungsten surface contribute to its fractional coverage 9, and 
this in turn leads to a work-function change A <p . But reference 11 bows to a tradition of 
thermionics and represents the arrival rate with a cesium reservoir temperature T Eg . 
When this reference uses a logarithmic term for effects of ion adsorption, the corre- 
lation fails at low coverages. Thus, the present work relies on the simple effective 
version (eqs. (10) and (17) of ref. 11): 
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A cp = 


-8.75 9 


1 + 1.77 9 


3/2 


0<|A<p| <2. 88V, 0 < 0 <0.6 


(5) 



/ 





1.37 - 1.34/ A<f> N 

\ 2 - Te 1- 

b + Tcs -< 

A + In 

7.02X10 _6 (1 - 9) 

1 

\8.75 e t 

1 * cs r 

11 606 11 606 1 

1 

H 

03 

to 

l 

1 - 9 


( 6 ) 


Between 9 values of 0.6 and 0. 8, A cp equals about 2. 88 volts. 

In reference 11 the constants A and B appear in an expression that correlates 
Taylor and Langmuir's vapor-pressure data for cesium in millimeters of mercury (mul- 
tiply by 133. 3224 for N/m 2 ): 


hi p cs = A - 



(7) 


The values of the constants A and B are given in the following table for various cesium 
reservoir temperatures. 


Cesium reservoir 
temperature, 

T Cs- 

K 

Constants in eq D (7) 

A 

B, 

K 

224 to 302 

17.70 

9300 

302 to 408 

16.17 

8820 

408 to 551 

15.76 

8660 

551 to 745 

15.38 

8440 


Only T Cg remains undefined; it depends on the arrival rate of equation (2): 

t C s= WM 2 

2 ™Cs« \ * / 


( 8 ) 


With these equations and the solution for the positive-ion sheath, the calculations 
form a standard iterative pattern: 

(1) Select temperatures for the electrode and plasma and assign the plasma electron 
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concentration. 

(2) Compute sheath results for <p = 4. 59 volts (refs. 1 to 4). 

(3) Determine the total cesium arrival rate (eq. (2)). 

(4) Calculate T Cg (eqs. (7) and (8)). 

(5) Obtain 9 and A cp (eqs. (5) and (6)). 

(6) Add A (p to 4. 59 volts to give a new <p. 

(7) Average the new and previous values of (p . 

(8) Compute sheath results for the new average cp. 

(9) Iterate steps (3) to (8) until <p changes by less than 0. 001 volt. 

This cycle with the separate sheath iteration in it produces complete emitter results and 
a detailed sheath description in about 0. 03 minute on an IBM 7090 computer. The com- 
plete calculation procedure appears in reference 16. 


DISCUSSION OF RESULTS 

Figures 1 to 5 evolve the effects of collisionless positive- ion sheaths on plane poly- 
crystalline tungsten surfaces in cesium plasmas. In parts (a) of the figures equilibrium 
results appear; for the (b) parts and figure 5 the plasma electron temperature runs up to 

200 K hotter than that of the rest of the system. Plasma concentrations derive from 

12 13 14 

thermal ionization by 10 ,10 , and 10 electrons per cubic centimeter. Tempera- 

tures range from 1700 to 2600 K giving pressures between 10“^ and 10 torr (multiply by 
133. 3 for N/m ). And the electrodes with bare work functions of 4. 59 volts emit ions, 
electrons, and atoms with effective work functions influenced by the sheath fields and 
adsorbed cesium on their surfaces. The results presented in figures 1 to 5 and any gen- 
eralization based on them apply only for these conditions with the restrictions given in 
the preceding sections and the previous publications on this sheath model (refs. 1 to 3). 
For the present results the selected temperatures prevent electrode melting, the Debye 
lengths are shorter than 1 mean free path for cesium charge exchange, and the net par- 
ticle current densities remain below 6 percent of their respective random plasma fluxes. 

Figure 1(a) shows overall sheath voltages for work functions of polycrystalline tung- 
sten with adsorbed cesium and for a constant work function of 4. 59 volts. If the work 
function were always 4. 59 volts, the sheath drop would rise with the plasma electron 
concentration at any fixed temperature. But actually the work functions climb or fall 
with increasing electron number densities depending on the temperature. For a given 
plasma electron concentration, the overall sheath voltage passes through a maximum as 
the system grows hot. With greater electron concentrations these extrema occur at high 
temperatures and appear lower and flatter. At high temperatures the sheath voltages 
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Figure 1. - Potential drops across positive-ion sheaths separating 
cesium plasmas (Saha) from plane polycrystalline tungsten sur- 
faces {(p = 4. 59 V). 
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System temperature, = T e p = Tjp = T a p, K 
(a) Equilibrium. 



Plasma electron temperature, T e p, K 
(b) With elevated plasma electron temperatures. 


Figure 2. - Cesium arrival rates at plane polycrystalline tungsten surfaces (<p - 4. 59 V) separated by 
positive-ion sheaths from plasmas (Saha). 


approach those for a 4. 59-volt work function. For these conditions, however, the true 
sheath drops never reach as high as 1 volt. Discussion of figures 2 to 5 indicates some 
reasons for this limit. 

Figure 2(a) reveals a similarity of atomic and total arrival rates only at low tem- 
peratures. When the system grows hotter, the arrivals of ions outnumber those of the 
atoms by over 10 to 1. Thus the cesium reservoir temperature generally fails to repre- 
sent cesium arrival rates properly in a thermionic diode operating at these conditions 
(ref. 10). 

Figure 3(a) gives fractional coverages by cesium of the polycrystalline tungsten sur- 
face. Initially the coverage drops off nearly exponentially with system temperatures; 
then it falls more and more rapidly with increasing temperature. Unfortunately only the 
high coverages here correspond to experimental data; the lower parts of these curves 
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(a) Equilibrium. (b) With elevated electron temperatures. 

Figure 3. - Fractional coverages of plane polycrystalline tungsten surfaces {<p = 4. 59 V) by cesium from 
plasmas (Saha). 

depend on extrapolations using the correlation of reference 11. 

Figure 4(a) indicates that work functions climb steeply with temperature initially, 
then level toward 4. 59 volts. Below the work functions lie the plasma electron potentials 
relative to the emitter Fermi levels. Reference 3 presents the equation for these equi- 
librium plasma potentials: 


+ 1.624J (9) 

The work function minus the plasma potential equals the overall sheath voltage plus 
the depression of the work function caused by the sheath field. This brings back the 
sheath-drop question of figure 1. Why do equilibrium sheath voltages for a constant 
work function drop steadily as temperatures rise while those for cesium on tungsten go 
through maxim ums? First, these sheaths build primarily with emission - the electron 
sheaths with emitted electrons, the positive- ion sheaths with emitted ions (refs. 1 to 3). 
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Electron voltage above Fermi level, V 



(a) Equilibrium. 



(b) With elevated electron temperatures. 

Figure 4. - Work functions for plane polycrystalline tungsten surfaces {<p = 4. 59 V) separated by positive-ion 
sheaths from cesium plasmas (Saha). 
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Second, ion emission grows greater with higher cesium arrival and work functions and 
with lower electrode temperatures (ref. 8): 


N iE 


N 


aE 


exp 

e(<P - I) 

_ * t e _ 


2 


( 10 ) 


But figures 2(a) and 4(a) show that small work functions hold ion production down while 
cesium arrivals are high and temperatures low. Then when work functions are high, 
low cesium arrivals and high temperatures drive the ion yield downward. As a result, 
sheath voltages pass through maximums where conditions are optimum for making and 
trapping ions. 

Figures 1(b), 2(b), 3(b), 4(b), and 5 present results for near-equilibrium positive- 
ion sheaths. Here the plasma electron temperature ranges up to 200 K higher than that 
of the rest of the system. Yet for parts (b) of the figures the general comments for the 
(a) parts still apply. Now however the temperature of the emitter, ions, and atoms re- 
mains constant at the initial value for each curve. The trends of parts (b) derive pri- 
marily from the decrease in cesium arrivals at the electrode as the plasma electron 
;emperature rises (fig. 2(b)). Because the plasma has a constant ionic concentration, 
ligher electron temperatures require fewer neutrals to maintain the fixed number of ions 
ref. 5); this causes the lower cesium arrival rate. 

As references 1 to 3 show, the plasma potentials for these positive-ion emission 
sheaths differ for equilibria and near equilibria. Equation (9) applies to equilibrium, and 
the following expression gives plasma potentials relative to emitter Fermi levels for 
near- equilibrium positive-ion sheaths: 
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Values from equation (11) appear in figure 4(b). 

Figure 5 also presents a distinct difference between the equilibrium and near- 
equilibrium positive-ion sheaths. Equilibria produce no net fluxes of particles, but net 
current densities climb abruptly from zero as the plasma electrons grow hotter than the 
rest of the system. After this sharp rise the curves seem to approach exponentials. 
Because the plasma electron temperature rises above that of the other particles and the 
electrode, these net currents flow from rather than to the plasma. 
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Figure 5. - Near-equilibrium currents through plane 
positive-ion sheaths separating polycrystalline tung- 
sten surfaces i(p = 4. 59 V) from cesium plasmas (Saha) 
with elevated electron temperatures. 

CONCLUDING REMARKS 

This study indicates effects of positive-ion sheaths on total cesium arrivals, surface 
coverages, and work functions of polycrystalline tungsten emitters. The results come 
from equilibria and systems with elevated plasma electron temperatures. Thus the 
present work parallels reference 1. That investigation, however, uses parametric work 
functions. In these calculations work functions vary below the 4. 59-volt value of the 
bare electrode depending on the amount of adsorbed cesium. This difference allows an 
evaluation of the influence of work functions on sheath structures. For comparison, data 
computed with a constant work function of 4. 59 volts also appear here. 

The results of the present study verify the strong effects of ionic arrivals on prop- 
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erties of electrodes with positive-ion sheaths in cesium plasmas (refs. 1, 2, 3, 10, 17, 
and 18). For cesium diodes the implications of this phenomenon are obvious. Positive- 
ion sheaths increase electron emission by reducing effective work functions with sheath 
fields and cesium adsorption. Proper selection of conditions can build positive-ion 
sheaths that produce more power. But understanding and controlling the internal mech- 
anisms of the diode depend on knowing the characteristics of the emitter. And a safe 
assumption at this point is that emitter properties have never been adequately determined 
for an efficiently operating thermionic diode (ref. 10). This state of ignorance persists 
both experimentally and theoretically. 

Lewis Research Center. 

National Aeronautics and Space Administration, 

Cleveland, Ohio, March 20, 1968, 

120-33-02-01-22. 
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conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof !’ 

— National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and 
technical information considered important, 
complete, and a lasting contribution to existing 
knowledge. 

TECHNICAL NOTES: Information less broad 
in scope but nevertheless of importance as a 
contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: 
Information receiving limited distribution 
because of preliminary data, security classifica- 
tion, or other reasons. 

CONTRACTOR REPORTS: Scientific and 
technical information generated under a NASA 
contract or grant and considered an important 
contribution to existing knowledge. 


TECHNICAL TRANSLATIONS: Information 
published in a foreign language considered 
to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information 
derived from or of value to NASA activities. 
Publications include conference proceedings, 
monographs, data compilations, handbooks, 
sourcebooks, and special bibliographies. 

TECHNOLOGY UTILIZATION 
PUBLICATIONS: Information on technology 
used by NASA that may be of particular 
interest in commercial and other non-aerospace 
applications. Publications include Tech Briefs, 
Technology Utilization Reports and Notes, 
and Technology Surveys. 


Details on the availability of these publications may be obtained from : 

SCIENTIFIC AND TECHNICAL INFORMATION DIVISION 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 



